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The well known effect of the linear damping on the moving nonlinear Schrodinger soliton (even 
when there is a supply of energy via the spatially homogeneous driving) is to quench its momentum 
to zero. Surprisingly, the zero momentum does not necessarily mean zero velocity. We show that two 
or more parametrically driven damped solitons can form a complex traveling with zero momentum 
at a nonzero constant speed. 

All traveling complexes we have found so far, turned out to be unstable. Thus, the parametric 
driving is capable of sustaining the uniform motion of damped solitons, but some additional agent 
is required to stabilize it. 



I. INTRODUCTION 



The amplitude of a nearly-harmonic wave propagating in a nonlinear dispersive medium satisfies a nonlinear 
Schrodinger equation. Confining ourselves to the generic, cubic, nonlinearity of the 'focusing' type, the resulting 
nonlinear Schrodinger equation is of the form 

i*t + Vxx + 2|*| 2 * = -17*; 7>0. (1) 

The —i'y'if! term in the right-hand side accounts for dissipative losses (which were assumed to be small in the derivation 
of eq.JTJ.) In the underlying physical system the dissipation is compensated by pumping the energy into the system, 
in one way or another. The pumping is modeled by adding a driving term to the right-hand side of eq.QJ. 

Like a simple pendulum, the distributed system can be driven externally or parametrically. The typical form of the 
corresponding amplitude equation is 

i^t + ^ xx + 2\^\ 2 ^ = he mt (2) 

and 

i#t + ^xx + 2|#| 2 $ = ti$e 2iat - Z7*, (3) 

respectively. (The overline in the right-hand side of © indicates complex conjugation.) Both the externally and 
parametrically driven nonlinear Schrodinger equations arise in a great variety of physical contexts. In particular, the 
parametric equation describes the nonlinear Faraday resonance in a vertically oscillating water tank jj]. Q and 
the effect of phase-sensitive amplifiers on solitons in optical fibers £| . The same equation controls the magnetization 
waves in an easy-plane ferromagnet exposed to a combination of a static and microwave field |5| and the amplitude 
of synchronized oscillations in vertically vibrated pendula lattices la. 

Both equations J2J an d © exhibit soliton solutions 0, IS0 , DT0 HI j stable and unstable 0, || , which can also 
form (stable and unstable) multisoliton complexes [Tot fill ~ll2fl. All localized solutions that have been found so far, 
were non-propagating. In fact, it is widely accepted that the nonlinear Schrodinger solitons simply cannot travel in 
the presence of dissipation. This perception is based on the rate equation 

P = -2 7 P, (4) 

which is straightforward from (J2J and J3J|. Here P is the total field momentum, 

if 00 — — 
p =~ ('fx^-^x^)dx. (5) 

* J -oo 

In the undamped case (7 = 0) the momentum is conserved; however if 7 > 0, P decays to zero and this seems to 
suggest that a solitary wave, initially moving with a nonzero velocity, will have to slow down and eventually stop |13| . 

Another indication that only quiescent solitons are possible in the damped-driven Schrodinger equation, comes 
ostensibly from the singular [2, Hj] and Inverse Scattering-based perturbation theory [3, Ha, Ha] ■ Here we should 
mention however that these techniques are well developed only in the one-soliton sector and in the case of several 
well separated solitons. They either make use of the smallness of the perturbation in the right-hand side of ©-© 
plPl Il5l| or utilize an explicit form of the perturbed soliton (to study its stability and bifurcation) [l^. In any case, 
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the resulting finite-dimensional system of equations for the parameters of the soliton and radiations, leads to the 
conclusion that the soliton's velocity has to decay to zero as t — > oo. 

Meanwhile, the moving solitary waves could play a significant role in a variety of physical situations modeled by 
the damped-driven nonlinear Schrodinger equations. Stable traveling waves could compete with non-propagating 
localized attractors; unstable ones might arise as long-lived transients and intermediate states in spatiotemporal 
chaotic regimes. Both types of moving solitary waves could mediate energy dissipation in damped-driven systems. 
One more reason for not rejecting the unstable solutions outright is their possible persistence within the (directly or 
parametrically driven) Ginzburg-Landau equations of which the Schrodinger equations ©-10 are special cases |18| . 
The diffusion and nonlinear damping (the terms ic\^ xx and — ic2\^\ 2n ^, to be added to the right-hand side of ©-©) 
are known to have a stabilizing effect on the Ginzburg-Landau pulses; hence the unstable Schrodinger solitons may 
gain stability as they are continued to nonzero positive c\ and c 2 . 

The purpose of this paper is to show that the damped-driven nonlinear Schrodinger equations do support solitary 
waves traveling with nonzero velocities. For the demonstration of this fact we confine our study to the 'parametrically 
driven Schrodinger only. The externally driven equation is left as an object of future research. 

Two complementary strategies will be pursued to achieve our goal. First, in section ITTll we consider the motionless 
damped solitons (V — 0, 7 7^ 0) and derive the condition under which they can be continued to nonzero velocity. 
Having identified values of 7 for which this condition is satisfied, we perform the numerical continuation obtaining 
a branch of solitary waves with nonzero V and 7. Our second approach is presented in section IfVI the idea is to 
continue undamped traveling waves (7 = 0, V ^ 0) to nonzero dampings. We show that this is only possible if the 
traveling wave has zero momentum. For complexes with P = 0, we then carry out the numerical continuation in 7. 
In section [3 we discuss the consistency of results obtained within these two complementary approaches. 

We examined, numerically, stability of all solutions obtained within both approaches. The general framework of the 
stability analysis is outlined in section [n] Results of this analysis are presented along with results of the numerical 
continuation. Finally, section Ivll summarizes conclusions of our study. 

II. MATHEMATICAL PRELIMINARIES 

For purposes of this paper we transform equation J2J to an autonomous form. First, we normalize the driving 
frequency O to unity; after that, the substitution #(x, t) — e lt ip(x,t) takes eq.(j2Jl to 

iipt + ipxx + Zlipl 2 ^ - iP = hip - ijip. (6) 

This is the representation of the parametrically driven damped nonlinear Schrodinger equation that we are going to 
work with in this paper. We confine ourselves to uniformly traveling solutions of the form 

Mx,t)=rl>(x-Vt)=1>($, (7) 

where "0(0 ~~ * as l£l — > 00 ■ These satisfy an ordinary differential equation 

-iVip s + + 2\ip\ 2 ip - if) = hn\> - iytp. (8) 

The analytical part of this paper deals mainly with identifying those of the previously found solutions of (JSJ 
with V = or 7 = which can be continued in V and 7, respectively. The actual continuation will be carried out 
numerically. Our numerical method employs a predictor-corrector continuation algorithm with a fourth-order accurate 
Newtonian solver. Typically, the infinite line was approximated by an interval (—100, 100). The discretization step 
size was typically 0.005. The numerical tolerance was set to be I0 -10 ; that is, the grid solution would be deemed 
accurate if the difference between the left- and right-hand sides in JSJl were smaller than I0~ 10 . 

Along with the continuation of solutions in V and 7, we will be analyzing their stability to small perturbations. To 
this end, we linearize equation eq.jfjj in the co-moving frame of reference. Letting ip(x,t) — u(£) + iv(£) + 5ip(£,t), 
where u and v are the real and imaginary part of the solution that we are linearizing about, and assuming that the 
linear perturbation depends on time exponentially: 

6j,(£,t) = e At [Su(0+iSv(0], 

we arrive at an eigenvalue problem 



ttoy= (A + 7)Jy, 



(9) 
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where the operator TLq is defined by 



. -<9f + 1 + h - 6u 2 - 2v 2 -Vd% - Auv . 

n °-\vdt-4uv -dl + l-h-6v 2 -2u 2 >> 



the skew-symmetric matrix J is 



-1 

1 



and the column vector y(£) = (Su, 5v) T . The eigenvalue problem @ was solved by expanding Su and Sv over a Fourier 
basis, typically with 500 modes, on the interval (—50,50). 

The last point that we need to touch upon in this preliminary section, is the integrals of motion of eq.|j2Jl, or, more 
precisely, the quantities which are conserved in the absence of dissipation. When 7 = 0, the equation @ conserves 
the momentum (given by eq.JHJl where one only needs to replace \f — > ip), and energy, 



E= / + \W ~ \W + h R.e^)eb- (11) 

J — oo 

In the damped case, the momentum decays according to the rate equation |J3J| while the energy satisfies 

E = 2j( f \ip\ 4 dx - E ] . (12) 



III. CONTINUATION OF DAMPED SOLITONS TO NONZERO VELOCITIES 



A. Continuability criterion 

Our first strategy is to attempt to continue stationary solutions with nonzero 7 to nonzero V. Two basic soliton 
solutions, denoted ip + and ?/>_, are available explicitly: 

ip± (x) = e~ l0± A± sech (A±x) ; (13) 
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A± = y/l± y /h*-'y 

1 7 7T 

6+ = -arcsin-, 6L = --0+. 

The two solitons can form a variety of stationary complexes. These are denoted, symbolically, tpi \, ipr_\ |_), 

ipi 1 ), and so on 12]. Let ^>o( x ) be a particular complex; we want to find out whether it can be continued in V. 

Assuming there is a solution ^>(£; V) such that ip(£; 0) = ^o(0 (= i'oix)), we expand V) in powers of V as 

m V) = e- w { uo(0 + iv (0 +V[ui(0 + ivi(0] + V 2 [u 2 (0 + iv 2 (0] + ...} , (14) 

where the constant phase will be chosen at a later stage. We also expand h and 7: h = ho+hiV+..., 7 = 70+71 V+.... 
Substituting into (JSJ, the order V 1 gives 

where the operator C has the form 



v i J \ ~ u o J \ v o 



r _ 1 -&x + 1 + ^0 cos 20 - 6uq - 2v 2 70 + h sin 20 - 4u t>o 1 ... 



-70 + fro sin 20 — 4w vo — ^ + 1 — ho cos 20 — 2uq — 6«q 

the constant matrix B is given by 



B 



— /iicos20 — 71 — /iisin20 
71 — /iisin20 /iicos20 
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and the primes over uq and Vq indicate derivatives with respect to x. (In i|15|) and i|16|l we have replaced £ with x as 
£ coincides with x for V = 0.) According to Fredholm's alternative, eci. (|15fl has a bounded solution ui(x),Vi(x) iff 
the vector in the right-hand side is orthogonal to the kernel of the Hermitean-conjugate operator £t; 

/ (2/ ' w) (-u' Q ) dx + j(y^)®( u v a ) dx = - (17) 

Here y(x) = (y, w) T is the eigenvector of & associated with the zero eigenvalue: C)y = 0. That the operator C) has a 
zero eigenvalue follows from the fact that the operator C has one — namely, the translation eigenvalue corresponding 
to the eigenvector (u' , v' ) T . Eq. l|17|l gives a necessary continuability condition of damped quiescent solitons to nonzero 
velocities. 



B. Noncontinuability of the 'building blocks' 

It is quite easy to check that when 70 7^ 0, the individual -0+ and tp^ solitons (the basic 'building blocks' of which 
all complexes are constructed) are not continuable to nonzero V. Choosing 9 = 9 + for the ip+ soliton and 9 = 9_ for 
the 0_ (where 9± are to be computed from the bottom formula in (| 1 3|) with 7 = 70 and h = /i ), we get Vq(x) = 0, 
70 — hosin29 = and so the 2x2 matrix £, ea. l|16|) . becomes upper triangular. The zero mode of & can now be 
readily found. 

Consider, for instance, the ip + case. The zero mode satisfies 



-dl + A\- Qui 
2 7o -d 2 x + A 2 _-2u 2 



0, 



hence y{x) = u' (x) and w(x) is found from 

{-dl + A 2 _ - 2u 2 )w = -2 1U ' (x). (18) 

Using the explicit expression for uo(x), uq(x) = A + sech(A + x), the operator in the left-hand side of l|18(l can be written 
as A 2 AL$ — e), where e = 2h$ cos{29+) / A 2 + ; Lq is given by 

L o = -d 2 x + 1 - 2sech 2 X, 

and X — A + x. The operator Lq has familiar spectral properties; in particular it has a single discrete eigenvalue 
Eq = associated with an even eigcnfunction zq = sechX, while its continuous spectrum occupies the semi-axis 
Ek > 1. Consequently, for < e < 1 (that is, for h < yl + To), the operator L — e is invertible and a bounded 
solution w(x) of l|18() exists and is unique. It can be found explicitly, but this is not really necessary for our purposes. 
All we need to know is that, since L is a parity-preserving operator, w(x) has the same parity as the right-hand side 
in l|18JI . i.e. it is an odd function. For that reason the second integral in equation (|17|l vanishes and the necessary 
continuability condition reduces to 



7 J u (x){L - e^u'^dx = 0. (19) 
This quadratic form can be easily evaluated by expanding u' (x) over eigenfunctions of the operator Lo: 

u' (x) = / U(k)z k (X)dk, 



where LoZk(X) = (1 + k 2 )zk(X). (The 'discrete' eigenfunction zo(^0 does not appear in the expansion as it has 
the opposite parity to u'q{x).) Utilizing the orthonormality of the eigenfunctions, the continuability condition i|19|) is 
transformed into 

As e < 1, this condition can obviously not be satisfied (unless 7 = 0). 

In the case of the ip- soliton the analysis is similar. In this case the continuability condition l|2(J|l is replaced by 



and this cannot be met for the same reason as ea. l|2U|l . 
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FIG. 1: A fragment of the bifurcation diagram for stationary multisoliton complexes (adapted from [13 ■) Shown is the energy 
mill of the complex as a function of h. The bottom branch pertains to symmetric two-soliton complexes ip(++) and ip(--) and a 

three-soliton solution ipr | \; the top branch includes the three-soliton states i/>( +++ ) and ipt ), as well as a five-soliton solution 

tpi | | y The thick curve corresponds to stable and thin curves to unstable solutions. The black dots indicate points where the 

integral H21H equals zero and therefore moving solitons are allowed to bifurcate off. 



C. Continuation of the complexes 



Turning to the complexes of the solitons ip+ and ip-, the phase of the complex varies with x and therefore the 
matrix C cannot be made triangular no matter how we choose the constant 9 in i|14[l . For this reason, aggravated 
by the fact that the multisoliton solutions are not available explicitly, the continuability condition i|17|l cannot be 
verified analytically. Resorting to the help of computer, we evaluated the eigenfunction y(x) associated with the zero 
eigenvalue of the operator & numerically. (Here we set 9 = 8 + ). 

All damped soliton complexes found in |l2j| . were symmetric; that is, the corresponding u and v are even functions 
of x. Therefore, the operator C) whose potential part is made up of u(x) and v(x), is parity preserving and all its 
eigenfunctions pertaining to non-rcpcatcd eigenvalues arc cither even or odd. As we move along a continuous branch 
of solutions, the parity of the eigenfunction has to change continuously. Since the parity equals either +1 (for even 
functions) or —1 (for odd functions), the only opportunity left to it by the continuity argument, is to remain constant 
on the entire branch. For that reason it is sufficient to determine the parity of the eigenfunction for one specific value 
of ft — and we will know it at all other points. Our numerical calculation shows that the eigenfunction y(x) is odd 
on all branches reported in |l2j. Consequently, the second term in (|17|) is always zero and we only need to evaluate 
the first term. 

The vanishing of the term involving coefficients fti and 71 in ea. (|17fl implies that it was not really necessary to 
expand h and 7 in powers of V. This fact has a simple geometric interpretation. As we will see below, for the fixed 
7 the continuable solutions occur only at isolated values of h; hence they exist only for ft and 7 lying on continuous 
curves in the (ft, 7)-plane. Each curve results from an intersection of some surface in the three-dimensional (ft, 7, V)- 
space and the V — 0-plane. The fact that one does not have to alter ft and 7 when continuing the solution to nonzero 
V, indicates that these surfaces are orthogonal to the V — plane along their curves of intersection. 

Having found the solution ip(x) — u(x) + iv(x) at representative points along each branch, we obtained the eigen- 
function y(x) at these points and evaluated what remains of the integral (|17f) : 

J(yv' -wu' )dx = I(h). (21) 

The integral / is a continuous function of ft, and it was not difficult to find points on the curve at which it changes 
from positive to negative values, or vice versa. 

We examined two branches of multisoliton solutions obtained previously ^3 (FigD- The integral 1(h) was found 
to change its sign at three points, marked by black dots in FigQ] (Although it may seem from the figure that / equals 
zero right at the turning points, in the actual fact zeros of / do not exactly coincide with the turning points.) We 
were indeed able to numerically continue our solutions in V from each of these three points. Results are presented in 
FigH (a)-(c). 

The point '1' in Fig ^corresponds to the stationary complex an d lies just above the turning point where the 

ip(++) turns into ij)r__y (The turning point has h = 0.83504217 while 1(h) = for ft = 0.8353.) This solution has four 
positive real eigenvalues in the spectrum of the associated linearized operator and hence is unstable. The "f(V) curve 
which results from the continuation of this solution in V, is shown in Fig[2Ia). As V grows from zero, the solution 
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FIG. 2: Bifurcation curves branching off the points marked by black dots in Fig0 The curves illustrate the relation between the 
value of the damping 7 and velocity V at which the wave may travel for that 7. Each curve begins at the point 7 = 0.565 on the 
vertical axis. The insets show representative solutions at internal points of each branch. (Solid line: real part; dashed line: imaginary 
part.) Note the logarithmic scale of V in (b). Here, and in all other diagrams, arrows indicate our direction of continuation. 



looses its even symmetry (see the inset to Fig|2Ja)) while the four positive eigenvalues collide, pairwise, and become 
two complex conjugate pairs with positive real parts. After reaching a maximum velocity of approximately 0.65, the 
curve turns back toward V = 0, with 7 first growing but then also turning toward 7 = 0. The solution transforms 

into a (strongly overlapped) ip(j ) complex. As V and 7 tend to their zero values, the separation between the ip+ 

and ip- constituent solitons in the complex grows to infinity. The spectrum becomes the union of the eigenvalues of 
the individual -0+ and tp- solitons; in particular, it includes a complex-conjugate pair with positive real part, and a 
positive real eigenvalue. Thus the entire branch shown in Fig|2|[a) is unstable. 

One more comment that we need to make here, concerns the validity of the continuation scenario presented in 
Figl2I a ) f° r other values of h. Note that if we chose a smaller value of 7 in Fig^ the value of h corresponding to the 
point '1' would also be smaller. (For example, for 7 = 0.548 the integral 1(h) vanishes at the point h = 0.82.) For 
this smaller h the final product of the continuation turns out to be not a pair of infinitely separated stationary ip+ 
and tp- but a totally different complex. This is discussed below in section [Vj see also FigsQJc) and[5j 

Another branch bifurcates off at the point marked '2' in Fig^ Here h = 0.867. The corresponding j(V) diagram 
is displayed in Fig^b). As we move along the branch departing from V = 0, the original stationary complex ipi \ 
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transforms into a solution displaying three widely separated peaks in its real part: one corresponding to a strongly 

overlapping complex ip( | y, the next one to the ip + and the last one to the tp_ soliton. After passing a turning 

point, the curve is reapproaching, tangentially, the V = 0-axis. However, having reached V = 2.2 x 10~ 8 , it suddenly 
turns back and the velocity starts to grow again. The separation between the solitons decreases and the solution can 

now be interpreted as a strongly overlapping four-soliton complex ^(++h ) (shown in the inset to Fig|5fb)). As we 

continue further, the four constituent solitons regroup into two complexes, ip(++) and ipr^ ). The distance between 

the two complexes grows rapidly and, for certain finite V and 7 (at the endpoint of the curve in FigEJb) ) becomes 

infinite. At this point we have two coexisting solutions, and 0(h )> an d so this point corresponds to the point 

of self-intersection of the curve shown in Fig|2fa). Continuing the two solutions, separately, from the endpoint of the 
curve in Fig|2fb), we reproduce the diagram of Fig^a) for a slightly different value of h (i.e. for h = 0.867.) 

The entire branch shown in Fig[3(b) is unstable. The start-off stationary solution ipi \ has three positive real 

eigenvalues in its spectrum; one of these persists for all V and 7 while the other two collide and form a complex- 
conjugate pair with a positive real part. 

The branch continuing from the point '3' in FigQ for which h — 0.863645, leads to the least expected solutions. 
The resulting j(V) curve is shown in Fig[2{c). For points lying on the 'spiral' part of the curve, the function ip(x) is 
equal to a constant in a relatively large but finite region, and zero outside that region. (See the inset to Fig|2fc).) The 
constant is ip^ = {A^ / \/2)e~ %e - ; it defines a stationary spatially uniform solution to eq.@. (This flat background 
is unstable with respect to the continuous spectrum perturbations. Figuratively speaking, our pulse solution ij)(x) 
represents a 'droplet' of the unstable phase in the stable one.) On one side (at the rear of the pulse) the zero 
background is connected to the background by a kink-like interface. In the front of the pulse, the interface has 

the character of a large-amplitude excitation, with the shape reminding the ip^ ) complex. As the curve j(V) spirals 

toward its 'focus' in Fig|5fc), the length of the region where ip{x) = is growing. The entire branch is unstable; the 

start-off ipi j solution already has two real positive eigenvalues in its spectrum and more appear as we move along 

the branch. Those additional positive eigenvalues are remnants of the 'unstable' interval of the continuous spectrum 
of the flat nonzero solution ip(°K 



IV. CONTINUATION OF TRAVELING WAVES TO NONZERO DAMPINGS 



A. Continuability conditions 

When 7 = 0, the equation JSJ has a plethora of localized solutions with nonzero V [l^, and our second strategy 
will be to attempt to continue these undamped traveling waves to nonzero 7. We start with establishing the necessary 
and sufficient conditions for such a continuation. 

A set of the necessary conditions can be easily derived using two integral characteristics of equation j6f . the 
momentum 

P = (i/2) J (V^ - M)dx, (22) 

and energy (|f ljl . No matter whether 7 equals zero or not, the uniformly traveling solitary waves (i.e. solutions of the 
form Q) satisfy P = E = 0. Using these relations in eqs.QJ and l|12|l with 7 7^ 0, we get 

P = 0, (23) 

and 

E = J \iP\ A dx. (24) 

Equations l|23|l - H24|) have to be satisfied by the undamped solutions continuable to nonzero 7. 

In fact, eas. <|23(l and l|24() are not independent. Indeed, multiplying Eq.(|SJ) by -0, adding its complex conjugate and 
integrating, gives an identity 

E — I l^dx = VP. (25) 



Letting P = in 125(1 . eq.J23J) immediately follows. Thus we can keep equation P = as the only necessary condition 
for the continuability to nonzero 7; eq.J23J is satisfied as soon as eq.(j22J is in place. 
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It turns out that P = is also a sufficient condition. To show this, we expand the field ip = u + iv in powers of 7: 

u = u + 7U1 + 7 2 u 2 + v = v + 7U1 + 7 2 «2 + — , 

substitute into JSJ and equate coefficients of like powers. (We could have also expanded h and V in 7, but, similarly 
to the continuation in V described in the previous section, the terms with coefficients hi and V\ cancel out of the 
resulting continuability condition.) At the order 0{^ 1 ), we obtain: 

Here the hermitian operator TIq is as in (|10|l where we only need to attach zero subscripts to u and v: 

_ ( -cf + 1 + h - Qui - 2v 2 -V0 6 - Au v \ 

Ho ~ \ Vdz - 4u v -d 2 + l-h-2u 2 -6v 2 ) ' {ll > 

Since the operator TCq has a zero eigenvalue, with the translation mode as an associated eigenvector, equation (j26(l is 
only solvable if its right-hand side is orthogonal to (u' , v' Q ): 

(u' ,v' )(~H^<% =0. (28) 

(Here the prime indicates the derivative with respect to £.) The expression in the left-hand side of l|28[) coincides with 
equation (1221) written in terms of the real and imaginary part of ip, and so the solvability condition l)28|) is simply 
P = 0. 

Now assume that P is equal to zero so that a bounded solution to equation (|26J) exists. All traveling waves found in 
|17| have even real and odd imaginary parts: uo(— x) = uo(x), vo(—x) = —vq(x). Noticing that the diagonal elements 
of the operator TIq are parity-preserving while the off-diagonal elements change their sign under the £ — > — £ reflection, 
we conclude that u±(x) is odd and vi{x) is even. 

Proceeding to the order C(7 2 ), we have 

n, ( u 2 \ _ / —v\ + u [6u 2 + 2v 2 } + 4u wii>i \ (0Q , 

The top entry in the right-hand side of l|29() is even and the bottom one odd; hence the right-hand side is orthogonal 
to the null vector (u' ,v' ) and a bounded solution U2{£,),V2(£,) exists. This time the u-component is even and the 
^-component odd: U2(— £) = u 2 (£), v z(~Q = ~ v ^(0- 

It is not difficult to verify that this parity alternation property guarantees the boundedness of u n (t;) and for 
all n. Therefore, equation JHJ has a localized solution — * as |£| — > 00) for sufficiently small 7. Thus if we have 
an undamped soliton traveling with zero momentum, it can be continued to nonzero values of 7. 



B. Continuable solutions: the bifurcation diagram of the undamped nonlinear Schrodinger 

In this subsection we review the P(V) dependence for the undamped solitons and solitonic complexes 01 • Of 
interest, of course, are points where the graph crosses the y-axis, i.e. where P(V) = 0. 

The simplest solutions arising for V — are, obviously, our stationary fundamental solitons and ?/>_ . These are 
given by eas. l|13|l where one only needs to set 7 = 0: 

ip+(x) — A + sech (A + x), ip-(x) = iA_scch (A-x), 

with A 2 ± = 1 ± h. Both ip + and ip- have zero momenta and therefore, are continuable to nonzero 7. However, the 
continuation does not produce any traveling waves in this case; all we get is our static damped solitons ip±, ea. (|13|l . 
Next, both ip + and ip~ admit continuation to nonzero V (for the fixed 7 = 0) 0|. As V is increased to 

c= y / 2 + 2\/r^ : ^, 

the width of the soliton rp_ increases, its amplitude decreases and the soliton gradually transforms into the trivial 
solution, tp = 0. On the resulting branch, the momentum vanishes only for V — and V — c and therefore, no 
damped branches can bifurcate off the traveling ip_ soliton. 



9 




FIG. 3: The momentum of the undamped traveling wave as a function of its velocity (a combined and advanced version of two 
diagrams from fl7|l. The dashed and dash-dotted curves pertain to the case of large driving strengths (here exemplified by h = 0.7). 
The starting point P = V = of the dashed curve corresponds to the stationary undamped V>+ soliton, which then transforms to the 
twist, then to a bound state of two twists and then to a complex of a twist and two tp- solitons. (This curve appeared in 17].) The 
dash-dotted offshoot is our new contribution to the diagram; it corresponds to an asymmetric solution, ipr_T), detaching from the 
VVtt) curve. The solid curve pertains to the case of small driving amplitudes (here h = 0.05). (This curve also appeared in [17|]) 
The points of its intersection with the P-axis correspond to stationary twist solitons; continuing each of these counterclockwise gives 
rise to a bound state of two tp+'s, while when continued clockwise each twist transforms into a complex of two ip-'s. More solution 
curves can be generated by the mapping V — > —V, P — > —P. 



We now turn to the soliton ip+. When h < 0.28, its fate is similar to that of the ip-: as V — > c, the soliton 
spreads out and merges with the zero solution. The momentum equals zero only at two points, V = and V = c; for 
< V < c, the momentum is positive. 

For h > 0.28, the transformation of the 0+ is more promising from the present viewpoint (see the dashed curve in 
FigEJ). As V is increased from zero, the momentum grows, then the branch turns back toward the V — axis. For 
some V < the momentum reaches its maximum and then decreases to zero. The point V = V\ where P(V\) = 
is of interest to us as a branch of damped solitons can bifurcate off at this point (and it really does, see subsection 
IIVDI ) Continuing beyond V%, the curve P(V) turns toward V = and then, after one more turning point, we have 
another zero crossing: P(V^) = 0. This is how far we have managed to advance in our previous work |l7j|. 

At this point we need to mention that the ip+ and ip- are not the only quiescent solitons for 7 = 0. The dashed 
P(V) curve in FigEl is seen to have one more intersection with the P axis, apart from the one at the origin. The 
corresponding solution represents a symmetric strongly overlapping complex of the if)+ and ip- solitons and was coined 
"twist" (symbolically ipx) in [l7|. The twist soliton arises both for h greater and smaller than 0.28. In the former 
region the twist obtains from the ^-continuation of the -0+ soliton while for h < 0.28, it is not connected to the 
(See the solid curve in FigEl) The continuation of the twist in V in the case h < 0.28 gives rise to a new branch of the 
undamped solutions which has a point of intersection with the P = axis, at some V — V\. A damped traveling wave 
is bifurcating off at this value of velocity; see the next subsection. We are using the same notation V\ in the small- 
and large-/i case in Fig0to emphasize the similarity of the resulting ~f(V) curves in the two cases (forthcoming). 

Returning to the case of large h, the entire dashed curve in Fig |3| corresponds to symmetric solutions: ip(— £) = 
i/ , (^). It turns out that there are also non-symmetric solutions; these were missed in |l7j . The real part of a non- 
symmetric solution is not even and imaginary part not odd. In particular, a pair of asymmetric solutions arise in a 
pitchfork bifurcation of the complex V'(tt) ; see the dash-dotted offshoot from the dashed curve in Fig|3| (The two 
asymmetric solutions are related by the transformation ip(£) — > £)> they obviously have equal momenta and hence 
are represented by the same curve.) Continuing the asymmetric branch we have the third zero crossing, at V = V3. 
When continued to positive P, the asymmetric solution acquires the form of a complex of ip- and i^t solitons, with 
the inter-soliton separation growing as P is increased. (Note that although the dashed and dash-dotted curve end at 
nearby points, they are not connected.) Our numerical analysis shows that branches of damped solitons do indeed 
detach at V\, V% and V3; these will be described in the next two subsections. 



C. Numerical continuation: Small driving amplitudes 

For small h, h < 0.28, our continuation departs from the twist soliton moving with the velocity V\ (the point of 
intersection of the solid curve with the horizontal axis in Fig|3|) The real part of this solution is even and imaginary 
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FIG. 4: Results of the numerical continuation of the undamped traveling solitons to nonzero 7. (a): small h; (b),(c): large h. The 
inset displays a representative solution at one of the internal points of the curve. (Solid line: real part; dashed line: imaginary part.) 
Each curve shown has a positive-velocity counterpart which arises by the mirror reflection V — » —V of the figure. 



part odd: ij)(—x) — ip(x). As we continue to nonzero 7, this symmetry is lost; a typical profile at the internal points 
looks like a non-symmetric complex of the ip- and and is displayed in the inset to FigQJa). The rest of Fig^Ja) 
shows the resulting "/(V) dependence. As 7 grows, the negative velocity of the traveling wave decreases in modulus. 
However the damping cannot be increased beyond a certain limit value; as we reach it, the 7(V r )-curve turns down 
(FigEt a ))- As V and 7 tend to zero, the separation between the ifj- an d ip+ solitons in the complex grows without 
bounds. 

These transformations of the solution are reflected by the behavior of the linearized eigenvalues in the eigenvalue 
problem 0. At the point V = Vi, 7 = of the ~f(V) curve, the twist solution has a quadruplet of complex eigenvalues 
±A, ±A which dissociates into two pairs of complex-conjugate eigenvalues Ai,Ai and A2,A2 (with ReAi < and 
Re A2 > 0) as 7 deviates from zero. As we move toward the maximum of the curve, the imaginary parts of Ai and A2 
decrease and the four complex eigenvalues move onto the real axis. At the point of maximum one of the resulting two 
positive eigenvalues crosses to the negative real axis, but the other one persists all the way to V = —0 and 7 = +0. 
Therefore the spectrum of eigenvalues on the 'downhill' portion of the curve is a union of eigenvalues of the ip- and 
-0+ solitons. The conclusion of the eigenvalue analysis is that the traveling complex whose bifurcation diagram is 
exhibited in Fig^Ja), is unstable for all V and 7. 
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D. Numerical continuation: Large driving amplitudes 

For h > 0.28 we have three starting points with P = corresponding to two intersections of the dashed curve and 
one of the dash-dotted curve with the horizontal axis in Fig|3J 

The j(V) curve emanating out of the point V\ is the top, arc-shaped, curve in Fig^Jb). For V = V\ and 7 = the 
solution is symmetric and its shape reminds two strongly overlapping twists. The linearized spectrum includes two 
complex quadruplets. As 7 deviates from zero, the symmetry is lost and the solution starts looking like an asymmetric 
complex of two pulses. The two complex quadruplets become four complex-conjugate pairs of eigenvalues, two with 
positive and two with negative real parts. Two of these pairs (one with Re A > and one with Re A < 0) move on to 
the real axis. After that one positive real eigenvalue crosses to the negative semi-axis, while the complex pair with 
Re A > crosses into the Re A < half-plane but then returns to Re A > 0. As V,7 — » 0, the separation between the 
ip— and solitons comprising this complex increases, and eventually the two constituents diverge to infinities. On 
the 'downhill' portion of the curve, the spectrum is a union of the spectra of the individual "0- and ip+ solitons; in 
particular, it includes a positive real eigenvalue and a complex quadruplet. Since there are eigenvalues with Re A > 
for all V, the entire branch is unstable. 

The second undamped traveling wave with zero momentum (point V2 on the bifurcation diagram Fig|3J) corresponds 
to a symmetric [ip(—x) — VK 2 -)] complex of two ip— and one twist soliton, symbolically ip(_ T _y The spectrum includes 
three complex quadruplets. As we continue in 7 and V , the symmetry is lost but the solution still looks like a complex 
of three solitons, see the inset to Fig0Jb). The bottom, spike-shaped, curve in FigQJb) depicts the corresponding ^{V) 
dependence. Unlike the branch starting at the value V — V\, this solution cannot be continued to zero velocities. 
Instead, the "f(V) curve turns back and, as 7 approaches zero from above, V tends to a negative value V4, with 
\Vi\ > I V2 1 - For sufficiently small 7 the corresponding solution consists of two ip- solitons and a twist in between, 
with the inter soliton separations growing to infinity as 7 — > 0, V — > V4. The associated eigenvalues perform rather 
complicated movements on the complex plane; skipping the details it suffices to mention that 'unstable' eigenvalues 
(real positive or complex with positive real parts) are present for all V. Hence the entire branch is unstable. 

Finally, the point V3 on the diagram FigOHrepresents two nonequivalent asymmetric solutions with zero momentum, 
^>i(£) and ^(C), with foiQ = ipi{— £)■ Consequently, there are two distinct 7(U)-branches coming out of this point 
(Figfflc))- One of these corresponds to a complex of two solitons; when continued to V = 0, it gives rise to the 
symmetric complex "0(++) with nonzero 7. (See the top curve in Fig^c)). Continuing the other asymmetric solution 
to V = 0, the corresponding value of 7 reaches a maximum at V ~ 0.3 and then tends to zero. (The bottom curve in 

FigEJc)). For sufficiently small V and 7 this solution represents a complex ipt (shown in the inset to FigQJc)). 

As V, 7 — > 0, the inter-soliton separation tends to infinity. Turning to the eigenvalues, the start-off solution at the 
point V3 has two complex quadruplets and a real positive eigenvalue in its spectrum. When we continue along the 
top curve in Fig^Jc), two complex eigenvalues move on to the positive real axis, so we end up with three positive 
eigenvalues. When we continue along the bottom curve, the movements of the eigenvalues are more involved but 
some of them always remain in the 'unstable' half-plane, Re A > 0. The upshot of the eigenvalue analysis is that both 
curves represent only unstable solutions. 

V. CONSISTENCY OF THE TWO APPROACHES 

To complete our classification of damped traveling solitons, we need to comment on what may seem to be an 
inconsistency between results obtained within the above two complementary approaches. The solution representing 
the well-separated ip+ and ip- solitons reported in sections II I II and IIVI can be reached by continuing both off the 
(7 = 0)- and (V — 0)-axes. (This branch connecting to the origin on the (V, 7)-plane appears both in Figs. Ufa) 
and^b).) Although such a curve should obviously not depend on the starting point of the continuation, one notices 

that the ip^ ) branches 'flowing into the origin' in Figs. EJa) and Efb) behave differently when traced backward 

(i.e. away from V = 7 = 0). While the curve in Fig[3[a) intersects the 7-axis, its counterpart in FigQJb) crosses 

the other, V-, axis. (Here the reader should not be confused by the fact that the i/jr^ ) branch in Fig^a) is shown 

for positive and its counterpart in Fig0Jb) for negative values of V. In view of the £ — > — £, V — > — V invariance 
of equation iJHJ , to each 7 there correspond two traveling waves, one with positive and the other one with negative 
value of V. Therefore, one should mirror-reflect Fig^Jb) prior to comparing it to FigJ^a). This reflection maps the 
solution of Fig^Jb) to the ip(+-) of FigJ^a).) 

To resolve the paradox, one needs to note that the two figures correspond to different values of h, FigJ^fa) to 
h = 0.8353 and FigQJb) to h — 0.7. It turns out that a qualitative change of behavior occurs for h somewhere 
between these values, more precisely between 0.82 and 0.8275: For h — 0.82 and smaller (in particular, for h = 0.7) 
the 7(V) curve has the form of an arc shown in Fig^Jb) (i.e. it crosses the V-axis as V is increased) while for 
h = 0.8275 and greater, the curve is already loop-shaped and does not reach to 7 = 0. This change of behavior, 



12 



0.6 

Y 

0.4 
0.2 


-0.2 
-0.4 

-0.6 — 

0.2 0.4 0.6 v 0.8 

FIG. 5: The comparison of the bifurcation diagrams for h — 0.82 and h — 0.8275. Solid curve: h — 0.8275; dashed curve: h — 0.82. 
Although the points of intersection of the dashed and solid curves with the V-ax\s do not coincide, i.e. 14(0.8275) 7^ 14(0.82) and 
Vi(0.8275) 7^ Vi(0.82), they are quite close to each other and hence we mark them as a single point (V3 and Vi, respectively.) The 
branches coming out of the point Vi on the V-axis are omitted for visual clarity. Note the shape of the dashed and solid curves near 
V ~ 0.7, 7 ~ 0.3, characteristic of phase portraits of 2D dynamical systems in the neighborhood of saddle points. 

accounting for the above 'inconsistency', is illustrated by Fig[S] which compares the j(V) dependences for h = 0.82 
and h = 0.8275. Fig[S]also serves to illustrate the different outcomes of the continuation of the complex ip(++) for 
h = 0.8353 and smaller h. (We note that for h = 0.8353 the continuation of the motionless vpr++) produces a pair 
of infinitely separated solitons ip + and ip- [Figl^a)] while for h — 0.7, the curve departing from the same type of 
starting point [i.e. from ?/'(++)] ends up at the undamped asymmetric solution traveling with nonzero velocity V3 
[FigHc)].) 

The above differences in behavior result from the presence of a saddle point on the (7, y)-plane, in the gap between 
the two lobes of the solid curve in Fig|SJ Indeed, the dashed and solid curves can be seen as sections of the surface 
h = h(j, V) by the horizontal planes h = 0.82 and h = 0.8275, respectively. The gap in the upper solid curve is then 
accounted for by letting h = ho + x 2 — y 2 in the vicinity of the gap. Here the constant ho lies somewhere between 
0.82 and 0.8275, and {x,y) is a pair of suitably chosen coordinates on the (7, y)-planc. 

VI. CONCLUSIONS 

One of the conclusions of this work is that by grouping into complexes, solitons (or, equivalently, solitary pulses) 
can adjust their total momentum to zero. By doing so they can travel with nonzero speed in the presence of damping 
— without violating the momentum decay law, P = — jP. Two identical solitons traveling at the same speed in the 
same direction have equal momenta; therefore, in order to arrange for P = the traveling complex inevitably has to 
include solitons of different varieties (i.e. both V'+' s an d tp-'s.) Consequently, the real and imaginary parts of the 
traveling complex will always be represented by asymmetric functions of £ = x — Vt. 

Although the possibility of non-decelerated motion may be out of line with the common perception of the soliton 
dynamics in weakly damped Hamiltonian equations, moving pulses are not unknown in strongly dissipative systems. A 
suitable example is given by the complex Ginzburg-Landau equation. Asymmetric Ginzburg-Landau pulses, uniformly 
traveling with nonzero velocities, were reported in [l9l ]. 

All moving solutions that we have found in this paper, turned out to be unstable. This instability admits a simple 
qualitative explanation — at least, for small dampings. In the undamped situation, the ip- solitons are unstable when 
traveling with small velocities while the ip+'s become unstable when moving sufficiently fast 01- ^ n the presence of 
dissipation the traveling wave has to include solitons of both varieties; on the other hand, the eigenvalues corresponding 
to small nonzero 7 should remain close to their (7 = 0)-counterparts. Therefore the spectrum of the traveling complex 
will 'inherit' unstable eigenvalues of either ^_ (for small velocities) or of the ip+ (for large velocities). 

Thus, despite the fact that the parametric driver can sustain the uniform motion of a damped soliton, an additional 
agent (such as, possibly, the diffusion and/or a nonlinear damping term) is required to make this motion stable. Here 
it is appropriate to refer, again, to the complex Ginzburg-Landau equation. Stable Ginzburg-Landau pulses arise 
as a result of a delicate balance of the whole series of terms, includ ing dispersion, cubic and quintic nonlinearity, 
diffusion, cubic gain and linear and quintic nonlinear damping [Tgt I20L l2l| . In a similar way, the gain/loss and 
spreading/steepening balances of the dampcd-driven traveling solitons could be restored by adding one or several 
missing agents. 
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